Topological Insulators at Room Temperature 
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Topological insulators are new states of quantum matter with surface states protected by the 
time-reversal symmetry. In this work, we perform first-principle electronic structure calculations for 
Sb2Te3, Sb2Se3, Bi2Te3 and Bi2Se3 crystals. Our calculations predict that Sb2Te3, Bi2Te3 and 
Bi2Se3 are topological insulators, while Sb2Se3 is not. In particular, Bi2Se3 has a topologically 
non-trivial energy gap of O.SeV, suitable for room temperature applications. We present a simple 
and unified continuum model which captures the salient topological features of this class of materials. 
These topological insulators have robust surface states consisting of a single Dirac cone at the F 
point. 

PACS numbers: 



Recently, the subject of time reversal invariant topo- 
logical insulators has attracted great attention in con- 
densed matter physics^^^^^'^'^'^^^^^^^^^^. Topological states 
of quantum matter are defined and described by the cor- 
responding topological terms in quantum field theory. 
For example, the quantum Hall effect is described by the 
topological Chern-Simons term^^. On the other hand, 
the electromagnetic response of three dimensional (3D) 
topological insulators are described by the topological 
term of the form Se = J d^xdtE ■ B, where E 

and B are the conventional electromagnetic fields and 
a is the fine structure constant. For a periodic system, 
all physical quantities are invariant under the shift of 
the parameter 6 by integer multiples of 27r. Therefore, 
all time reversal insulators, including strongly correlated 
and disordered ones, fall into two disconnected classes, 
described^ either by 6> = or by 6> = tt. Topological in- 
sulators are defined by ^ = tt, and this term physically 
described the topological magneto-electric effect^. Topo- 
logical insulators have surface or edge states with an odd 
numbers of gapless Dirac points. 

The 2D topological insulator displaying the quantum 
spin Hall (QSH) effect was first predicted for the HgTe 
quantum wells^. Recently, the edge state transport has 
been experimentally observed in this system^. The elec- 
tronic states of the 2D HgTe quantum wells are well de- 
scribed by a 2 + 1 dimensional Dirac equation where the 
mass term is continuously tunable by the thickness of 
the quantum well. Beyond a critical thickness, the Dirac 
mass term of the 2D quantum well changes sign from 
being positive to negative, and a pair of gapless helical 
edge states appear inside the bulk energy gap. This mi- 
croscopic mechanism for obtaining topological insulators 
by inverting the bulk Dirac gap spectrum can also be gen- 
eralized to other 2D and 3D systems. The guiding prin- 
ciple is to search for insulators where the conduction and 
the valence bands have the opposite parity, and a "band 



inversion" occurs when the strength of some parameter, 
say the spin-orbit coupling, is tuned. For systems with 
inversion symmetry, a method based on the parity eigen- 
values of band states at time reversal invariant points 
can be applied^. Based on this analysis, the Bi^Shi-x 
alloy has been predicted to be a topological insulator for 
a small range of x, and recently, surface states with an 
odd number of crossings at the fermi energy has been 
observed in angle-resolved photo-emission spectroscopy 
(ARPES) experimentsia. 

Since BixShi-x is an alloy with random substitutional 
disorder, its electronic structures and dispersion relations 
are only defined within the mean field, or the coherent po- 
tential approximation (CPA). Its surface states are also 
extremely complex, with as many as five or possibly more 
dispersion branches, which are not easily describable by 
simple theoretical models. Alloys also tend to have im- 
purity bands inside the nominal bulk energy gap, which 
could overlap with the surface states. Given the impor- 
tance of topological insulators as new states of quantum 
matter, it is important to search for material systems 
which are stoimetric crystals with well defined electronic 
structures, preferably with simple surface states, and 
describable by simple theoretical models. In this work 
we focus on layered, stoimetric crystals Sb2Tes^ Sb2Ses, 
Bi2Tes and Bi2Se^. Our theoretical calculations predict 
that Sb2Tes^ Bi2Tes and Bi2Se^ are topological insula- 
tors while Sb2Ses is not. Most importantly, our theory 
predicts that Bi2Ses has a topologically non-trivial en- 
ergy gap of O.SeV, therefore, it is a topological insu- 
lator at room temperature. The topological surface 
states for these crystals are extremely simple, described 
by a single gapless Dirac cone at the k = F point. We 
also propose a simple and unified continuum model which 
capture the salient topological features of this class of ma- 
terials. In this precise sense, this class of 3D topological 
insulators share the great simplicity of the 2D topological 
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FIG. 1: Crystal Structure (a) Crystal structure of Bi2Se3 
with three primitive lattice vectors denoted as ti,2,3- A quin- 
tuple layer with Sel — Bil — Se2 — Bil' — SeV is indicated 
by the red box. (b) Top view along the z direction. Trian- 
gle lattice in one quintuple layer has three different positions, 
denoted as A, B and C. (c) Side view of the quintuple layer 
structure. Along z direction, the stacking order of Se and 
Bi atomic layers is the sequence • • • — C{SeV) — A[Sel) — 

B{Bil) - C{Se2) - A{Bir) - BiSel') - C{Sel) . Sel 

(Bil) layer can be related to Sel' (Bil') layer by an inver- 
sion operation of which Se2 atoms play the role of inversion 
center. 



insulators realized in the HgTe quantum wells. 

Band structure and parity analysis. Bi2Se^^ Bi2Tes^ 
Sb2Tes, and Sb2Ses share the same rhombohedral crys- 
tal structure with the space group D^^ {R3m) with five 
atoms in one unit cell. We take Bi2Ses as an example 
and show its crystal structure in Fig. la, which has lay- 
ered structures with triangle lattice within one layer. It 
has a trigonal axis (three fold rotation symmetry), de- 
fined as z axis, a binary axis (two fold rotation symme- 
try), defined as x axis, and a bisectrix axis (in the re- 
flection plane), defined as y axis. The material consists 
of five-atom layers arranged along z direction, known as 
quintuple layers. Each quintuple layer consists of five 
atoms with two equivalent Se atoms (denoted as 5'el 
and Sel' in Fig. lb), two equivalent Bi atoms (denoted 
as Bil and Bil' in Fig. lb), and a third Se atom (de- 
noted as Se2 in Fig. lb). The coupling is strong be- 
tween two atomic layers within one quintuple layer but 
much weaker, predominantly of the van der Waals type, 
between two quintuple layers. The primitive lattice vec- 
tors 1 1^2,3 and rhombohedral unit cells are shown in Fig. 
1(a). Se2 site plays the role of inversion center and un- 
der inversion operation, Bil is changed to Bil' and 5'el 
is changed to 5'el'. The existence of inversion symmetry 
enable us to construct eigenstates with definite parity for 
this system. 

Ab initio calculations for £'62^63, Sb2Se^^ Bi2Tes and 



FIG. 2: band structure, Brillouin zone and parity 
eigenvalues. Band structure for Bi2Se3 without spin-orbit 
coupling(SOC) (a) and with SOC (b). The dashed line in- 
dicates Fermi level, (c) BZ for Bi2Se3 with space group 
R3m. The four inequivalent time-reversal invariant points 
are r(0, 0, 0),L(7r, 0, 0), F(7r,7r,0) and Z(7r,7r,7r). The blue 
hexagon shows the 2D BZ of projected (1,1,1) surface, in 
which the high-symmetry k points F, K and M are labeled, 
(d) The parity of the band at F point for the four materials 
Sb2Te3, Sb2Se3, Bi2Se3 and Bi2Te3. Here we show the par- 
ities of fourteen occupied bands, including five s bands and 
nine p bands, and the lowest unoccupied band. The product 
of the parities for the fourteen occupied bands is given in the 
bracket on the right of each row. 



Bi2Ses are carried out in the framework of PBE-type^ 
generalized gradient approximation(GGA) of the density 
functional theory (DFT)^^'^'^. ESTATE package^^ with 
plane- wave pseudo-potential method is used with k-point 
grid taken aslOxlOxlO and the kinetic energy cut- 
off fixed to 340eV. For Sb2Te3, ^^2^63 and ^^2^63, the 
lattice constants are chosen from experiments^^, while for 
Sb2Ses, the lattice parameters are optimized in the self- 
consistent calculation for rhombohedral crystal structure 
(a = 4.076A, c = 29.830A), due to the lack of experiment 
data. 

Our results are consistent with the previous 
calculations^^' In particular, we note that Bi2Ses 
has an energy gap about 0.3eV, which agrees well with 
the experimental data (about 0.2 — 0.3eV)^^'^^. In the 
following, we take the band structure of Bi2Ses as an 
example. Fig. 2 (a) and (b) show the band structure 
of Bi2Ses without spin-orbit coupling (SOC) and with 
SOC, respectively. By comparing the two figures one 
can see clearly that the only qualitative change induced 
by turning on SOC is an anti-crossing feature around 
F point, which thus indicates an inversion between the 
conduction band and valence band due to SOC effect, 
suggesting Bi2Ses to be a topological insulator. To 
firmly establish the topological nature of this material, 
we follow the method proposed by Fu and Kane^ and 



calculate the product of the parities of the Bloch wave- 
function for the occupied bands at all the time-reversal 
invariant momenta r,F, L,Z in Brillouin zone (BZ). 
As expected, we find that at T point the parity of one 
occupied band is changed upon turning on SOC, while 
the parity remains unchanged for all occupied bands 
at other momenta F, L, Z. Since the system without 
SOC is guaranteed to be a trivial insulator, we conclude 
that Bi2Ses is a strong topological insulator. The same 
calculation is performed for the other three materials, 
from which we find that Sb2Tes and Bi2Tes are also 
strong topological insulators, and Sb2Ses is a trivial 
insulator. The parity eigenvalues of the highest 14 bands 
below the fermi level and the first conduction band at 
r point are listed in Fig. 2 (d). From this table we can 
see that the product of parities of occupied bands at F 
point changes from the trivial material Sb2Ses to the 
three non- trivial materials, due to an exchange of the 
highest occupied state and the lowest unoccupied state. 
This agrees with our earlier analysis that an inversion 
between the conduction band and valence band occurs 
at F point. 

To get a better understanding of the inversion and the 
parity exchange, we start from the atomic energy levels 
and consider the effect of crystal field splitting and spin- 
orbit coupling to the energy eigenvalues at F point, which 
is summarized schematically in three stages (I), (II) and 
(III) in Fig. 3 (a). Since the states near Fermi surface are 
mainly coming from p orbit als, we will neglect the effect 
of s orbitals and starting from the atomic p orbitals of 
Bi (65^6j9^) and Se (4s^4p^). In stage (I), we consider 
the chemical bonding between Bi and Se atoms within 
a quintuple layer, which is the largest energy scale in the 
current problem. First we can recombine the orbitals in 
a single unit cell according to their parity, which results 
in three states (two odd one even) from each Se p orbital 
and two states (one odd one even) from each Bi p orbital. 
The formation of chemical bonding hybridize the states 
on Bi and Se atoms, thus push down all the Se states 
and lift up all the Bi states. In Fig. 3 (a), these five 
hybridized states are labeled as |P1^^^^^), |P2^^^^^) and 
|P0~^^), where the superscripts +, — stand for the par- 
ity of the corresponding states. In stage (II), we consider 
the effect of the crystal field splitting between different 
p orbitals. According to the point group symmetry, the 
Pz orbital is split from px and py orbitals while the latter 
two remain degenerate. After this splitting, the energy 
levels closest to the Fermi energy turn out to be the Pz 
levels I PI J) and |P2~). In the last stage (III), we take 
into account the effect of SOC. The atomic SOC Hamil- 
tonian is given by Hgo = XI ■ S^ with l^S the orbital 
and spin angular momentum, and A the SOC parameter. 
The SOC Hamiltonian mixes spin and orbital angular 
momenta while preserving the total angular momentum, 
which thus leads to a level repulsion between |P1+, |) and 
|P1^+^^, i), and similar combinations. Consequently, the 
|P1+, I (I)) state is pushed down by the SOC effect and 
the |P2~,| (I)) state is pushed up. If the SOC is large 
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FIG. 3: Band Sequence, (a) Schematic picture of the 
evolution from the atomic Px,y,z orbitals of Bi and Se into 
the conduction and valence bands of Bi2Sez at F point. 
The three different stages (I), (II) and (III) represent the ef- 
fect of turning on chemical bonding, crystal field splitting 
and spin-orbit coupling, respectively (see text). The blue 
dashed line represents the Fermi energy, (b) The energy lev- 
els |PlJ) and |P2~) of Bi2Se3 at F point versus an artifi- 
cially rescaled atomic spin-orbit coupling X(Bi) = xXo(Bi) — 
1.25xeV, A(5'e) = xAo(5'e) = 0.22xeV (see text). A level 
crossing occurs between these two states at x = Xc — 0.6. 



enough (A > Ac), the order of these two levels is reversed. 
To see this inversion process explicitly, we also calculate 
the energy levels |PlJ) and |P2~) for a model Hamilto- 
nian of Bi2Ses with artificially rescaled atomic SOC pa- 
rameters X{Bi) = xXo{Bi)^ X{Se) = xXo{Se)^ as shown in 
Fig. 3 (b). Here Xo{Bi) = 1.25eV and Ao(5'e) = 0.22eV 
are the realistic value of Bi and Se atomic SOC param- 
eters, respectively^^ From Fig. 3 (b) one can see clearly 
that a level crossing occurs between |PlJ) and |P2~) 
when the SOC is about 60% of the realistic value. Since 
these two levels have opposite parity, the inversion be- 
tween them drives the system into a topological insulator 
phase. Therefore, the mechanism for the 3D topological 
insulator in this system is exactly analogous to the mech- 
anism in the 2D topological insulator of HgTe. In sum- 
mary, through the analysis above we find that Bi2Ses 
is topologically nontrivial due to the inversion between 
two Pz orbitals with opposite parity at F point. Similar 
analysis can be carried out on the other three materials, 
from which we see that Sb2Tes and Bi2Tes are qualita- 
tively the same as Bi2Se^^ while the SOC of Sb2Tes is 
not strong enough to induce such an inversion. 

Topological surface states. The existence of topologi- 
cal surface states is one of the most important proper- 
ties of the topological insulators. To see the topological 
features of the four systems explicitly, we calculate the 
surface states of these four systems based on ah initio 
calculation. First we construct the maximally localized 
Wanneir function(MLWF) from the ah initio calculation 
using the method developed by N. Marzari et al?'^-^. 
With these MLWF hoping parameters, we employ itera- 
tive method^^^^^ to obtain the surface Green function of 
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FIG. 4: Surface states. Energy and momentum depen- 
dence of the local density of states (LDOS) for (a) Sb2Se3, 
(b) Sb2Te3, (c) Bi2Se3 and (d) Bi2Te3 on the [111] surface. 
Here warmer color represents higher LDOS. The red regions 
indicate bulk energy bands and the blue regions indicate bulk 
energy gap. The surface states can be clearly seen around 
r point as red lines dispersing in the bulk gap for Sb2Te3, 
Bi2Se3 and Bi2Te3- No surface state exists for Sb2Se3. 



the semi- infinite system. The imaginary part of the sur- 
face Green function is the local density of states (LDOS), 
from which we can obtain the dispersion of the surface 
states. When calculating the surface Green function, we 
only use the bulk's MLWF hoping parameters as semi- 
infinite systems without considering surface corrections. 
Due to the layered structure of these materials, we ex- 
pect the surface reconstruction effect to be minor for 
[111] surface. The surface LDOS on the [111] surface 
for all the four systems are shown in Fig. 4. For £'62^63, 
Bi2Ses and Bi2Tes^ one can clearly see the topologi- 
cal surface states which form a single Dirac cone at T 
point. On comparison Sb2Ses has no surface state and 
is a topological trivial insulator. Thus the surface state 
calculation agrees well with the bulk parity analysis, and 
confirm conclusively the topologically nontrivial nature 
of the three materials. For Bi2Ses the Fermi velocity of 
the topological surface states is vf — 5.0 x lO^m/s, which 
is similar to that of the other two materials. 

Low energy effective model. Since the topological na- 
ture is determined by the physics near F point, it is 
possible to write down an simple effective Hamiltonian 
to characterize the low-energy long-wavelength proper- 
ties of the system. Starting from four low lying states 
|P1+,T (i)) and |P2-,t (j)) at F point, such a Hamil- 
tonian can be constructed by the theory of invariants'^ 
for the finite wavevector k. Based on the symmetries of 
the system, the generic form of the 4x4 effective Hamil- 
tonian can be written down up to the order of O(k'), 
and the tunable parameters in the Hamiltonian can be 
obtained by fitting the band structure of our ab initio 
calculation. The important symmetries of the system 
are time-reversal symmetry T, inversion symmetry / and 



three fold rotation symmetry Cs along the z axis. In the 
basis of (|P1+, T) , \P2-, T> , i> , |P2ja)), the rep- 

resentation of the symmetry operations are given by T = 
/C-zcr^ 012x2, I = l2x2 0T3 and Cs = exp (^f cr^ (8) 12x2), 
where JC is the complex conjugation operator, a^'^'^ and 
^x,y,z (^ejiote the Pauli matrices in the spin and orbital 
space, respectively. By requiring these three symmetries 
and keeping only the terms up to quadratic order in k, we 
obtain the following generic form of the effective Hamil- 
tonian: 



H{k) = eo(k)l4x4 
/Mik) A,k, 
■ ■ -M{k) 
A2k+ 






A2k+ 




o(k') (1) 



with k-\ 



± iky, eo(k) = C + Dik^ + D2kj^ and 



A^(k) = M — Bikl - B2k\. By fitting the energy spec- 
trum of the effective Hamiltonian with that of the ab 
initio calculation, the parameters in the effective model 
can be determined. For Bi2Ses, our fitting leads to 
M = 0.28eV, Ai = 2.2eV-A, A2 = 4AeV-A, Bi = lOeV- 
A', B2 = 56.6eF-A2, C = -0.0068eF, Di = l.SeV- A^, 
D2 = 19.6eF-A^. Except for the identity term eo(k), the 
Hamiltonian ([1]) is nothing but the 3D Dirac model with 
uniaxial anisotropy along z direction and k dependent 
mass terms. From the fact M,Bi,B2 > we can see 
that the order of the bands |T1+, T (i)) and |T2-, j (j)) 
are inverted around k = compared with large k, which 
correctly characterizes the topologically non-trivial na- 
ture of the system. Such an effective Dirac model can be 
used for further theoretical study of the Bi2Ses system, 
as long as the low energy properties are concerned. For 
example, as one of the most important low energy proper- 
ties of the topological insulators, the topological surface 
states can be obtained from diagonalizing the effective 
Hamiltonian ([1]) with an open boundary condition, with 
the same method used in the study of two-dimensional 
quantum spin Hall insualtor— . For a surface perpendic- 
ular to the z direction {i.e., [Ill] direction), the surface 
states are described by a 2 x 2 massless Dirac Hamiltonian 



Hsurf^kx 1 ky^ 





A2k^ 



A2k. 





(2) 



in the basis of (|k, |) , |k, |)). Here the surface 
state wavefunction |k, |(|))is a superposition of the 
|Pl+,T(i)) and |P2+,T(i)), respectively. For A2 = 
4.1eV-Aobtained from the fitting, the fermi velocity of 
surface states is given by = ^2/^ — 6.2 x lO^m/s, 
which agrees reasonably with the ab initio results shown 
in Fig. 4 (c). In summary, the surface effective theory ([2]) 
characterizes the key features of the topological surface 
states, and can be used in future to study the surface 
state properties of the Bi2Ses family of topological insu- 
lators. 

In conclusion we have theoretically predicted a new 
class of topological insulators Sb2Tes, Bi2Tes and 
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Bi2Ses- In particular, Bi2Ses has a large topologically 
non-trivial energy gap ~ 0.3eV, sufficient for room tem- 
perature operation. The topologically nontrivial nature 
of these three materials originates from a band inversion 
at r point, similar to the strained 3D HgTe^'^^ and two- 
dimensional HgTe quantum wells^. The topologically ro- 
bust surface states are studied by ab initio method, which 
consist of a single Dirac cone around the T point. We 
have also obtained a 4 x 4 effective theory to characterize 
the bulk properties at low-energy and long- wavelength, 
and a 2 X 2 massless Dirac model to describe the surface 
states. 

The topological surface states can be directly veri- 
fied by various experimental techniques, such as ARPES 
and scanning tunneling microscopy (STM). In the re- 
cent years, evidences of surface states have been observed 
for Bi2Se3 and ^^2^63 in ARPES^^ and STM^^ exper- 
iments. In particular, the surface states of Bi2Tes ob- 
served in Ref^ had a similar dispersion as we obtained 



in Fig. 4 (d), which were also shown to be quite stable 
and robust, regardless of photon exposure and temper- 
ature. Thus this experimental result strongly supports 
that the surface states have topological origin. Further 
experimental studies on the surface state properties, such 
as ARPES, STM and transport measurements are nec- 
essary to verify our prediction. Moreover, the 3D topo- 
logical insulators are predicted to display the universal 
topological magneto-electric effecti^ when the surface is 
coated with a thin magnetic film. Compared with the 
Bii-xSbx alloy, the surface states of the Bi2Ses family of 
topological insulators contain only a single fermi pocket, 
making it easier to open up a gap on the surface by mag- 
netization and to observe the topological Faraday/Kerr 
rotation^ and image magnetic monopole effect^^. If ob- 
served, such effects can be unambiguously identified as 
the experimental signature of the non-trivial topology of 
the electronic properties. 
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